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Lifshitz and Excited State Quantum Phase Transitions in Microwave Dirac Billiards 
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We present experimental results for the density of states (DOS) of a superconducting microwave 
Dirac billiard. Such systems serve as models for the electronic properties of finite graphene sheets. 
The DOS exhibits two sharp peaks which evolve into van Hove singularities with increasing system 
size. These divide the Fermi surface into regions governed by the relativistic Dirac equation and 
by the non-relativistic Schrodinger equation, respectively. We demonstrate that there a topological 
transition appears as a neck-disrupting Lifshitz transition in the number susceptibility and as an 
excited state transition in the electronic excitations. Furthermore, we recover the finite-size scaling 
typical for excited state quantum phase transitions involving logarithmic divergences. 
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Introduction.- Graphene, a monolayer of carbon 
atoms forming a hexagonal lattice, has attracted a lot of 
attention in recent years due to its extraordinary prop- 
erties associated with the shapes of the conduction and 
the valence band. These touch each other conically, thus 
implying a linear dispersion relation. As a consequence, 
close to these so-called Dirac points excitations are de- 
cribed by a Dirac Hamiltonian [l| . The inset of Fig. [1] 
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FIG. 1: (Color on line) Density of states p(/) obtained 
from the resonance spectrum (see text below) of a microwave 
Dirac billiard (black line). The red line results from a tight- 
binding model [3]. The region of low density around the 
Dirac frequency fo ~ 23.36 GHz is delimitated by two sharp 
peaks, van Hove singularities, at f~jj = 21.98 GHz and 
f^u — 24.87 GHz. The inset shows isofrequency lines of the 
numerically determined Fermi surface of propagating modes 
in the quasi-momentum plane {qx,qy) (see text). 



shows isofrequency lines [Si] of the band structure (black 
Hues), i.e., the Fermi surface in the plane of the quasi- 
momentum vector components (qx,qy). The Dirac (K) 
points are located at the corners of the Brillouin zone 
(BZ). In their vicinity the isofrequency lines form circles, 
that deform into triangles with the distance from them. 



This relativistic regime is bordered by saddle points at 
the M points. At the centre of the Brillouin zone, the 
F point, the conduction (valence) band has a maximum 
(minimum) . In its vicinity the isofrequency lines form cir- 
cles and the Fermi surface has a parabolic shape. There 
the non-relativistic Schrodinger equation applies. Ac- 
cordingly, the Fermi surface can be separated into two 
independent relativistic regions and one non-relativistic 
one. The change of the topology takes place at the M 
points, where the Dirac cones merge into the paraboli- 
cally shaped surface. At the energies of the M points 
the density of states p (DOS) diverges logarithmically in 
an infinitely extended graphene sheet Q. These are the 
so-called van Hove singularities, which were predicted to 
exist in general in 2-dimensional crystals with a periodic 
structure §■ 

The topological phase transition associated with these 
singularities may be realized either by changing the lat- 
tice structure [61, '7l| or by applying a chemical potential, 
see Ref. [S] ■ We will demonstrate in this letter that it can 
be identified with a neck-disrupting ground-state Lifshitz 
transition !|| . Furthermore we will show that it can also 
be interpreted as an excited state quantum phase transi- 
tion (ESQPT) in the single-particle excitations. Transi- 
tions of the latter type have been observed in the equiv- 
alent bosonic system. A particularly close analogy with 
the present case is provided by the 2-dimensional vibron 
model [lOi] describing transverse vibrations of molecules, 
but ESQPTs have been studied in many other models 
as well 11 1|- 14| , including a detailed study of quasi-spin 
systems [1^. We investigated properties of the DOS of 
graphene experimentally in its microwave analog and re- 
covered this finite-scaling behavior, thus corroborating 
its universal nature. Figure [T] presents the density of 
states (DOS) of the microwave Dirac billiard. Its exper- 
imental determination is described in the next section. 

Experiment.- Microwave billiards have been used 
since two decades as analog systems for the study of non- 
relativistic quantum phenomena |16l . Il7l [. The present 



letter is about experiments which were performed with 
microwave photonic crystals. Generally, photonic crys- 
tals [la, lly| are the optical analog of a solid and the 
frequencies of wave propagation as function of the two 
components of the quasi momentum exhibit a band struc- 
ture. The microwave photonic crystals were composed of 
metallic cylinders constituting a triangular lattice |20l.[21| 
and squeezed between two metal plates. The structure of 
the crystal's first two bands is similar to the Fermi sur- 
face of graphene, that is, it is Dirac like in the vicinity 
of their touch points [22|. There, the scalar Helmholtz 
equation describing the sytem reduces to the Dirac equa- 
tion of massless spin-1/2 particles. The linear dispersion 
relation in graphene is due to the presence of two inter- 
penetrating triangular lattices with threefold rotational 
symmetry in the hexagonal lattice .23 .] . In the triangu- 
lar photonic crystal the honeycomb lattice is formed by 
the voids enclosed by, respectively, three cylinders. We 
use this analogy to experimentally investigate relativistic 
phenomena occurring in graphene. Various effects have 
already been studied as,_e.g., pseudo-diffusive trans por t 
near the Dirac point jM-iM] , the quantum Hall effect [23] , 
Zitterbewegung [IJ, and edge states [H, [H, ill . 

Here we present results associated with the properties 
of the DOS determined experimentally for two supercon- 
ducting Dirac billiards '26!]. They were constructed by 
enclosing a photonic crystal in a rectangular microwave 
billiard. The rectangular billiard basins containing the 
metallic cylinders were milled out of a brass plate with 
side lengths 420.0 x 249.4 mm'^. One Dirac billiard con- 
tained 267 metallic cylinders and had the lattice constant 
aL = 20mm, the other one 888 with a^ = 12mm. The 
radius of the cylinders was R = aL/4. Figure [2] dis- 
plays a photograph of the Dirac billiard with 888 cylin- 
ders. The brass lids and the basin were lead coated to 




FIG. 2; (Color on line) Photograph of the superconducting 
microwave Dirac billiard containing 888 metal cylinders. It 
is constructed from brass and coated with lead. The lid is 
shifted with respect to the billiard body. 

achieve superconductivity at liquid helium temperature. 
To ensure a good electrical contact the lids were screwed 
tightly to each cylinder. The height of the Dirac billiards 
was h — 3mm. Hence, up to the maximum frequency 
of 50 GHz, only the lowest transverse magnetic mode 
with the electric field vector perpendicular to the top 
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FIG. 3: A measured transmission spectrum of the microwave 
Dirac billiard depicted in Fig. [2] (upper panel) . It is termi- 
nated by two stop bands, where no wave propagation is pos- 
sible. The lower panel shows a zoom into the region of par- 
ticularly low resonance density around the Dirac frequency. 



and bottom plates was excited. Accordingly, the vecto- 
rial Helmholtz equation reduces to a scalar one which 
is mathematically identical to the Schrodinger equation 
of the corresponding 2-dimensional multiple-scattering 
problem. 

Microwave power was coupled into and out of the res- 
onator via wire antennas, that reached a few millimeters 
into the resonator through holes in the lid, for the mea- 
surement of the resonance spectra. A Vector Network 
Analyzer (VNA) measured the relative phase and ampli- 
tude of the output to the input signal. In the upper panel 
of Fig. |3] a transmission spectrum of the Dirac billiard 
with 888 cylinders measured in the frequency region be- 
tween 19.5 GHz and 30.5 GHz is depicted. It is bordered 
by two stop bands corresponding to the band gaps in the 
band structure where no wave propagation is possible. 
Furthermore, we observe a region with an exceptionally 
low resonance density around the Dirac frequency of the 
Dirac points. The lower panel of Fig. |3] shows a zoom into 
that region. Transmission spectra were measured with all 
possible combinations of two out of a total of five anten- 
nas attached to the lid at different positions. Since the 
resonances had high-quality factors Q > 5 • 10^, we could 
resolve all resonances and determined 1651 eigenfrequen- 
cies. 

Figure [T] shows the DOS p{f) obtained from the mea- 
sured resonance spectra, i.e., the number of states per 
frequency interval versus the frequency of the centre of 
the interval. We chose a frequency interval of 100 MHz. 
The red curve was computed based on a tight-binding ap- 
proach [2] which takes into account not only the nearest- 



neighbor, but in addition the second-nearest neighbor 
coupling and also the overlap integral between the wave 
functions located at neighboring voids [2a] ■ In the region 
of low density the DOS vanishes linearly around the Dirac 
frequency fo = 23.36 GHz with \f — fo] ^ 0- It corre- 
sponds to the relativistic region in the Fermi surface of 
graphene, i.e., the propagation of electromagnetic waves 
is governed by the Dirac equation [22]. This region is 
delimitated by two sharp peaks at /^ — 21.98 GHz and 
/J}^ = 24.87 GHz, the van Hove singularities i5]. In the 
frequency range below f~^ and above /^ the Helmholtz 
equation reduces to the Schrodinger equation of the cor- 
responding quantum billiard [la, [l3|. This defines the 
non-relativistic region. A closer look at the experimental 
DOS reveals, that the amplitudes and the typical fre- 
quencies of the oscillations of the experimental DOS are 
smaller in the frequency range between the two van Hove 
singularities than below and above. This already indi- 
cates that both regions are governed by different wave 
equations. 

At the van Hove singularities the DOS diverges loga- 
rithmically. This happens strictly only in an infinitely ex- 
tended graphene sheet or photonic crystal. In the Dirac 
billiards used in the experiments, however, the sharp 
peaks at /^ have a finite height p"^^^. We determined 
it for the experimental DOS of the two microwave Dirac 
billiards, and also performed numerical studies with pho- 
tonic crystals of various sizes and shapes. For a compari- 
son of these results we rescaled the frequencies such that 
the distance /J}^ — f~fj between the van Hove singulari- 
ties, i.e., the group velocity, was the same for all systems. 
We chose the rescaling / — 7> / such that /J^ — /J^ = 2. 
The experimental and numerical studies revealed that 
the maxima of the DOS, p'^^^, or rather those of the 

renormalizcd DOS, n"^^^ = 



- &-J^ -max behave like 



aNc (In(iVc) + b) 



(1) 



with Nc the number of unit cells, i.e., of hexagons formed 
by the voids in the photonic crystal. The quantities a and 
b are fit parameters, where the latter depends on the ge- 
ometry of the billiard, i.e., is system dependent, while 
the former takes a similar value a ^ 0.145 — 0.155 for all 
cases, i.e., seems to be universal. This finite-size scaling, 
typical of an ESQPT [l^ , and the fate of the isofrequency 
lines at the saddle points (see inset in Fig. [T|) suggests a 
description in terms of a neck-disrupting Lifshitz tran- 
sition [9]. We note, however, that all properties of the 
DOS that we observe in photonic crystals can also be ob- 
tained from the DOS for the vibrations perpendicular to 
the plane of an hexagonal lattice, as shown by Hobson 
and Nierenberg [23] in 1953. Except for the choice of the 
frequency scale their result coincides with the analytical 
expression for the DOS ^4;] investigated in the following 
two sections. 

Neck-disrupting Lifshitz transition in graphene.- In 
order to illustrate the relation between the van Hove 
singularity in the DOS and the neck-disrupting Lif- 



shitz transition in the analogous fermionic band struc- 
ture of graphene, we computed the number susceptibil- 
ity from the particle-hole polarization function, the Lind- 
hard function [30]. For this we used the simplest tight- 
binding model, which takes into account only nearest- 
neighbor hopping of strength i [2, |j] ■ Many aspects con- 
cerning the electronic excitations in graphene at weak 
coupling [j] can be studied analytically with this model 
to exemplify more general effects [34]. The retarded 
polarization function U^ (ui , p; fj,) , with w the excitation 
frequency, p the vector of momentum transfer and p, 
the chemical potential p33l . |34| . is a sum of particle- 
hole transitions within the same band, H:J, i.e., intra- 
band transitions and those arising from interband tran- 
sitions between the two bands, H;??. The polarization 
function yields for the (zero-temperature) number sus- 
ceptibility [33 



X — lim lim H 
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Hence, it coincides with the DOS p{lS) per area A of the 
graphene sheet at the Fermi surface to = p. It should be 
noted that only intraband transitions contribute. 

In Ref. [4| the zero of the DOS, identified with the 
Dirac point, is located at /x = 0, the van Hove singulari- 
ties are a,t p — zLt and the band gaps start at p — ±3i. 
When the chemical potential is chosen near one of the van 
Hove singularities we readil;j^ obtain from the analytical 
expression Eq. (14) in Ref. 



3JVc / 
^ ~ 2Tr^ At V 



lnlzl + 21n2 + C'(z; 



(3) 



where we introduced the parameter z = [\p\ — t)/t. Note 
the logarithmic divergence for z — > 0, which may be re- 
lated to the neck-disrupting Lifshitz transition in two di- 
mensions [3, l35| . 

As in Dirac billiards, the susceptibility or DOS does 
not diverge in a graphene sheet of finite area. To see how 
the heights of its maxima scale with the area we used 
periodic boundary conditions and integrated Eq. ([3|) over 
a small interval Az = {2tt)'^/Nc around the singularity. 
We again rescaled the energies such that the distance 
between the maxima equals two. Then we obtain for the 
height of the maxima of the renormalized DOS 
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-^ iV, In iV, - 2 In TT + 1 + ©(l/A^e) 



Note that ^ 



(4) 
0.15, thus confirming the experimental 
and the numerical findings, c.f., Eq. ([Ij. Thus the height 
of the maxima of the susceptibility at the van Hove singu- 
larities scales as tx"^'^^ = n"^^^/A ^ IniVc, in accordance 
with the finite-size scaling of a neck-disrupting Lifshitz 
transition. 

The divergence of x as z — >■ is caused by the infinite 
degeneracy of ground states of the 2-dimensional system 
when the Fermi surface passes through a van Hove sin- 
gularity. In the thermodynamic sense this can be con- 
sidered as a zero-temperature quantum phase transition 



with control parameter |/z|. However, unlike the cases of 
first- or second-order phase transitions, the susceptibility 
does not diverge with a power law in z but logarithmi- 
cally. Furthermore, the transition is due to a change of 
topology of the Fermi surface with no order parameter in 
the strict sense. 

Excited state quantum phase transition in the elec- 
tronic excitations.- In the following we will present a 
quasi-order parameter and demonstrate that, as in the 
equivalent bosonic system [10l - ll5l | , the singularity of the 
single-particle DOS as function of the excitation fre- 
quency can also be interpreted as an ESQPT in the 
spectrum of particle-hole excitations. For this we ana- 
lyzed the polarization function at zero-momentum trans- 
fer, Il^{uj,p — 0; /i). The associated spectral distribution 
Pph(a;) is given by {36| 



Pph('^) = ^-1 fim —^lmll^{uj,p;n) 



2ttp 



with the normalization 



Z — lim 



cLu - — — Im II''^ {ijj,p;pL) 
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(5) 



(6) 



It can be separated into contributions Z^ from intra- 
band, respectively, Z- from interband transitions. Here, 
Z+ is related to the /-sum rule [3J, |33, [Sgl , which ex- 
presses the conservation of the number of particles in 
terms of an identity. Indeed, Z^ is fixed due to charge 
conservation via this rule in terms of the 2-dimensional 
charge carrier density Hc and mass m, Z^ = ^. Near 
the centre of the Brillouin zone we have ric = p|,/(47r) 
and m ~ \/3Nc/ {tA) and the contribution from inter- 
band transitions behaves as Z-{ii) ~ TMIir^"^^ ~ \lAY ■ 
Hence it is suppressed such that Z{p) sa Z+(/i), and the 
sum rule is readily verified. 
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This approximation holds as long as one restricts to intra- 
band transitions and to the non-relativistic Fermi liquid 
regime either below or above the two van Hove singulari- 
ties, i.e., for |/i| > t. Near the Dirac cone, where |/i| <C i, 
on the other hand, the intraband transitions yield 



Z+{p) 




(8) 



with n'^ — /i^/(27rw|,). Thus, there the intraband /-sum 
rule scales with the square root of the carrier density n'^ 
relative to half filling. However, the contribution of the 
interband transitions to Z , 



Z^(m) 



■nt 
24^3 



ImI 



(9) 



can no longer be neglected for |/i| ^ t. Note, that the 
sum Z of the contributions Eqs. ([S]) and ©.is indepen- 
dent of pL and hence of the carrier density [3J, [33] . 



From these observations we conclude that the /-sum 
rule or Z[pL) can serve as a quasi-order parameter for 
the Lifshitz transition, indicating relativistic behavior 
for |/i|/i < 1 with Z K, const., as compared to the 
non-relativistic Fermi-liquid regime for |/i| > i, where 
Z(/i) « ^+(/i) decreases almost linearly with jjl. We ver- 
ified analytically that the derivative of Z[^i) with respect 
to pL diverges logarithmically at the Lifshitz transition, 
fj, = t. This reflects a singular behavior of the carrier 



density similar to that of x in Eq. ([3]), since Z^ 



for 




FIG. 4: (Color on line) Spectral distribution pph(aj) of verti- 
cal particle-hole excitations derived from the imaginary part 
of the polarization function as function of the rescaled chem- 
ical potential p/t and excitation frequency u/t showing the 
ESQPT at cj/i = 2 for \p\/t < 1. Also displayed in the front 
panel for a; = is the number susceptibilty x ~ pCa*)/^ to 
indicate the corresponding ground-state Lifshitz transitions 
at /i = ±i. 

We obtained the full spectral distribution pp\^{Lu) of 
particle-hole excitations Eq. ([5]) from explicit analyti- 
cal expressions for the polarization function. The result 
for H:^(a;, p; n) from intraband transitions was first de- 
rived in |33l |. Here we are particularly interested in the 
imaginary part of Il^{uj,p; n) for the interband transi- 
tions which reflects the ESQPT in the Dirac-cone re- 
gion for |/x| < t. For small momenta we computed the 
spectral density as described in Ref. ^g]. The result is 
summarized in Fig. 2] where we display pph(w). In the 
front panel we have included the number susceptibility 
x(/i) — p{fi)/A, to indicate the two ground-state Lifshitz 
transitions at fi — ±t. As in the Dirac cone approxi- 
mation, interband contributions to /9ph(w) vanish when 
w < 2/i because the vertical particle-hole excitations are 
then Pauli blocked. The ESQPT in the particle-hole exci- 
tation spectrum at zero momentum transfer is clearly vis- 
ible in Fig. m at w = 2/i. There the spectral distribution 
exhibits a logarithmic divergence which is directly linked 



to that of the single-particle DOS at the van Hove singu- 
larity Eq. ([2). Below this, for uj <2t and fj, <t, we have 
the relativistic behavior of the low-frequency excitations 
as discussed above. Above the ESQPT, for to > 2t, the 
density of particle-hole excitations decreases fast with in- 
creasing frequency, as it does in the normal Fermi-liquid 
regime. 

Conclusions.- We determined the DOS in high- 
precision experiments with two superconducting Dirac 
billiards and recovered the finite-size scaling govern- 
ing such transitions experimentally and also numerically 
with Dirac billiards of different sizes and geometries. The 
DOS we obtain is similar to that of transverse vibrations 
of an hexagonal lattice and, most importantly, to that of 
the electronic band structure of finite graphene sheets. 



In the second part of the letter, we show that the prop- 
erties of the observed DOS can be described in terms of 
a QPT and ESQPT arising from the topological Lifshitz 
neck-disrupting phase transition. We plan now to in- 
vestigate the features of the oscillations observed in the 
experimental DOS in Fig. [TJ These obviously differ in 
the frequency range between the van Hove singularities, 
where the wave propagation is governed by the relativis- 
tic Dirac equation, and that, where the non-relativistic 
Schrodinger equation applies (see Fig.[T]). The differences 
in the amplitudes and the frequencies of the oscillations 
will be extracted from their Fourier transforms. 
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